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According the Collapsar model long gamma-ray bursts (LGRBs) involve relativistic jets that 
puncture the envelope of a collapsing star, and produced the y-rays after they break out. This 
model provides a theoretical framework for the well known association between LGRBs and mas- 
sive stars. However although this association is supported by a wealth of observations, to this date 
there is no direct observational evidence for the emergence of the jet from the star. In other words 
there is no direct evidence for the Collapsar model. Here we show that a distinct signature of the 
Collapsar model is the appearance of a plateau in the duration distribution of the prompt GRB 
emission at times much shorter than the typical breakout time of the jet. This plateau is evident 
in the data of all major GRB satellites, and provides a direct evidence supporting the Collapsar 
model. It also enables us to place limits on the sizes and masses of LGRB progenitors; suggests 
the existence of a large population of choked (failed) GRBs; and indicates that the 2 s duration 
commonly used to separate Collapsars and non-Collapsars holds for BATSE and possibly Fermi 
GBM GRBs, but it is inconsistent with the duration distributions of Swift GRBs. 
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1. Introduction 

The connection of long GRBs (LGRBs) to collapsing massive stars is strongly supported by 
a wealth of observational evidence. The most notable evidence is the association of half a dozens 
GRBs with spectroscopically confirmed broad-line Ic supernovae (SNe), and the identification of 
"red bumps" in the afterglows of about two dozens more, which shows a photometric evidence of 
underlying SNe. The model that provides the theoretical framework of this association is known 
as the Collapsar model (MacFadyen & Woosley, 1999; MacFadyen et al., 2001). According to 
this model, following the core collapse of a massive star, a bipolar jet is launched at the center of 
the star. The jet drills through the stellar envelope and breaks out of the surface before producing 
the observed y-rays. However, although this model is supported indirectly by the LGRB-SN asso- 
ciation, to this date we could not identify a clear direct observational imprint of the jet-envelope 
interaction, thus there is no direct confirmation yet of the Collapsar model. 

Here we show that under very general conditions the time that the jet spends drilling through 
the star leads to a plateau in the duration distribution of the GRBs at times much shorter than the 
breakout time of the jet. This plateau exists in the duration distribution of all major GRB satellites, 
and provides a strong observational support for the Collapsar model. An analysis of this plateau 
also (i) supports the hypothesis of compact stellar progenitors, (ii) implies the existence of a large 
population of chocked jets that fail to break out of their progenitor stars and (iii) enables us to 
determine statistically the fraction of non-Collapsars from the total GRB sample and the threshold 
duration that separates Collapsars from non-Collapsars. Specifically it shows that this time differs 
from one satellite to the other. 

2. The propagation of a jet in the stellar envelope 

The jet propagates in the star by pushing the stellar material in front of it, leading to the for- 
mation of a "head" of shocked matter at its front. The head propagates at sub relativistic velocities 
along most of the star (Matzner, 2003; Zhang et al., 2003; Morsony et al., 2007; Mizuta & Aloy, 
2009; Bromberg et al., 2011a), even though the jet is ejected at relativistic velocities. Therefore as 
long as the head is inside the star it dissipates most of the jet energy, and it needs to be constantly 
supported by the relativistic jet material to propagate. This implies that there is a minimal amount 
of time, tj,, that the jet engine needs to operate to get a successful jet breakout (Bromberg et al., 
2011b): 



where L iJ0 is the isotropic equivalent jet luminosity, 6 is the jet half opening angle and we have 
used typical values for a long GRB. and are the radius and the mass of the progenitor star, 
where we normalize their value according to the typical radius and mass inferred from observations 
of the few supernovae (SNe) associated with long GRBs. If the activity time of the engine t e < t\, 
the jet fails to escape and a regular long GRB is not observed. 

After the jet breaks out of the star it expands and produces the observed y-ray emission at large 
distances from the stellar surface. In most GRB models (e.g. Sari & Piran, 1997) the observed 
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duration of the GRB, ty, reflects the activity time of the central engine after the break out of the jet 

ty = t e —tf,. (2.2) 

The distribution of ty is therefore a convolution of the distributions of the engine activity time and 
the breakout time combined with cosmological redshift effects. 



3. The duration distribution of Collapsars 



Under very general conditions, Eq. results in a flat distribution of ty at durations sig- 
nificantly shorter than the typical breakout time. A simple way to show that is by considering a 
single value of % for all Collapsars and ignoring, for simplicity, cosmological redshifts and detector 
threshold effects. The probability that a GRB has a duration ty equals, in this case, to the probability 
that the engine operating time is ty+tb, i.e. 

Py(ty)dty = p e (t h +ty)dty, (3.1) 

where p r is the probability distribution of the observed durations and p e is the probability distri- 
bution of the engine operating times. It is clear that p e (tb + ty) ~ p e (tb) = const for any ty <C %■ 
Moreover, if p e (t e ) is a smooth function that does not vary rapidly in the vicinity of t e « t;„ over 
a duration of the scale of %, then the constant distribution is extended up to times ty < ty. In the 
case of interest % and t e are determined by different regions of the star: the breakout time is set 
by the density and radius of the stellar envelope at radii > 10 10 cm, while t e is determined by the 
stellar core properties at radii < 10 8 cm. The core and the envelope are weakly coupled (Crowther, 
2007) and the engine is unaware whether the jet has broken out or not. Therefore, it is reasonable to 
expect that p e (t e ) is smooth in the vicinity of t e tj, and Pyity) sa const for ty < %. In the opposite 
limit, where % <C ty, then ty ss t e . Together eq. 3T reads: 



I PeVy) ty > h 

Remarkably, as we show in Bromberg et al. (2011c), the plateau exists also in the more gen- 
eral case, when tj, varies and the redshift distribution and detector thresholds are considered. The 
constant value of %, in this case, should be replaced by a "typical" breakout time, %, modulo red- 
shift corrections. This plateau in the GRB duration distribution is independent of specific details 
of the relevant distributions and in particular of the details of p e . It is a direct prediction of the 



Collapsar model that follows immediately from Eq. |2.2l Given an average redshift (1 + z) ~ 3 
and Eq. |2.1| , we expect the effective (redshift corrected) typical breakout time to be of order 50 
s. At long observed durations, ty 3> %, the distribution p y is determined by a convolution of p e 
and the distributions of the breakout times and bursts' redshifts. However p r cannot be flatter than 
p e , otherwise p e would dominate the distribution. Therefore, an extrapolation of p r (ty S> %) to 
durations shorter than % provides a lower limit to the number of events with t e < %. Namely, it 
is an estimate of the minimal number of choked bursts. At very short durations (less than 2 s) an 
additional population of shorter and harder GRBs (SGRBs) appears (Kouveliotou et al, 1993). It 
is well established that most SGRBs are not associated with death of massive stars (Nakar, 2007, 
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Figure 1: The T 90 distribution, dN/dT 90 , of BATSE (red), Swift (blue) and Fermi GBM (green) GRBs. 
Also plotted is the distribution of the soft (HR < 2.6) BATSE bursts (magenta). For clarity the Swift values 
are divided by a factor of 5 and the Fermi GBM by 15. At shorter times the sample is dominated by non- 
Collapsars. Note that the quantity dN/dT is depicted and not dN/dlogT as traditionally shown in such plots 
(e.g., Kouveliotou et al., 1993). The black lines show the best fitted flat interval in each data set: 5 — 25 s 
(BATSE), 0.7 — 21 s {Swift), and 2.5 — 31 s (Fermi). The upper limits of this range indicate a typical breakout 
time of a few dozens seconds, in agreement with the prediction of the Collapsar model. The distribution at 
times > 100 s can be fitted by a power law with an index —4 < a < —3. The soft BATSE bursts show a 
considerably longer plateau (0.4 — 25 s), indicating that most of the soft short bursts are in fact Collapsars 
(Bromberg et al., 201 lc). 

and references there in), namely they are non-Collapsars, and the above argument of a flat distribu- 
tion doesn't apply to them. Therefore, when considering the overall burst duration distribution we 
expect a flat section for durations significantly shorter than 50 s down to the duration where these 
non-Collapsars dominate. 

4. The observed distribution of the prompt GRB durations 

The observed duration of a GRB is characterized by T 9 o « tJl +z), during which 90% of the 
fluence is accumulated. Fig. 1 depicts the observed distribution of T$q for the three major GRB 
satellites: BATSE 1 , Swift 2 and Fermi GBM 3 . Note that we plot the quantity p 7 (T 90 ) = dN/dT 
and not dN/dlogT traditionally shown in such plots (Kouveliotou et al., 1993). As predicted, a 
plateau around 2-30 s is clearly seen in all distributions. The extent of the plateau varies slightly 
from one detector to another. This is expected given the different detection threshold sensitivities 
in different energy windows (see below). The regions marked in solid bold lines in Fig. 1 are 
consistent with a flat distribution with a x 2 P er degree of freedom of 0.6, 1.3,0.7 for BATSE, Swift 
and Fermi respectively (see Bromberg et al., 201 lc, for details). At the high end of the plateau the 

'http://swift.gsfc.nasa.gov/docs/swift/archive/grb_table 

2 http://gammaray.msfc.nasa.gov/batse/grb/catalog/current/ 

3 http://lyra.berkeley.edu/grbox/grbox.php 
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Tgo distribution decreases rapidly and can be fitted at long durations (MOO s) by a power law with 
an index, a, in the range — 4 < a < —3. 

To test our hypothesis that the plateau in the duration distribution indeed reflects the distri- 
bution of Collapsars and is not, for example, a consequence of adding two separate distributions: 
one decreasing (SGRBs) and one increasing (LGRBs), we also plot the distribution, dN/dT, of soft 
BATSE bursts (magenta). Since non-Collapsars are harder on average (Kouveliotou et al., 1993), 
we expect this sample to be less contaminated and the plateau to extend to shorter durations than 
in the whole BATSE sample. Considering the soft bursts as bursts with hardness ratio 4 , HR < 2.6, 
the median value of bursts with Tgo > 5 s, we find that the plateau in this sample extends from 25 
s down to 0.4 s (1.3 # 2 /dof), an order of magnitude lower than the original 5 s in the complete 
sample. This lands a strong support to the conclusion that the observed flat distribution is indeed 
indicating on the Collapsar origin of the population. It also implies that HR is a good indicator that 
effectively filters out a large number of non-Collapsars from the BATSE GRB sample. 

5. The fraction of non-Collapsar SGRBs 

At short durations, the GRB distribution is dominated by the non-Collapsars SGRBs (Nakar, 
2007, and references there in). These bursts are hard to classify since all their hard energy properties 
largely overlap with those of the Collapsars (Nakar, 2007). The least overlap is in the duration 
distributions and hence, traditionally a burst is classified as a non-Collapsar if T$o < 2 s. Even 
though this criterion is based on the duration distribution of BATSE it is widely used for bursts 
detected by all satellites. 

The fraction of Collapsars at short durations can be estimated by extrapolating the plateau in 
the duration distribution to short times. Since there is also an overlap with the SGRBs at long 
durations, the real heights of the plateaus are somewhat lower than what is shown in fig. 1. In 
Bromberg et al. (201 Id) we make a joint fit of the non-Collapsars and the Collapsars by modeling 
the duration distribution of the non-Collapsars as well. The fitted distributions are then used to 
calculate the fraction / of the non-Collapsar from the total number of observed GRBs as a function 
of Tgo. This fraction represents the probability that an observed GRB with a given duration in not 
a Collapsar. Fig. 2 depicts / as function of Tgo for BATSE, Swift and Fermi GBM samples. Table 
1 gives T$o at several selected values of this fraction. It can clearly be seen that using Tgo < 2 s 
as a method to identify non-Collapsars works reasonably well for BATSE and is marginal in the 
case of Fermi, but in Swift data it results in a large over estimate of the number of non-Collapsars. 
Adopting, for example, / = 0.5 as the threshold fraction that separates the two populations, we find 
that for BATSE the transition occurs at ~ 3.5 s, for Swift it occurs at ~ 0.6 s and for Fermi at ~ 1.7 
s (Note, however, that the errors in the Fermi data are quite large due to lack of statistics). This 
shift in the transition time is expected since non-Collapsars are on average harder than Collapsars 
and different detectors have different energy detection windows. BATSE has the hardest detection 
window, making it relatively more sensitive to non-Collapsar GRBs. Swift has the softest detection 
window making it relatively more sensitive to Collapsar GRBs. 

4 HR is the fluence ratio between BATSE channels 3 (100-300 keV) and channel 2 (50-100 keV). 
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Figure 2: The fraction / of non-Collapsars from the total number of observed GRBs as a function of the 
observed duration time, 790, in the (from top to bottom) BATSE, Swift & Fermi GBM samples. The shaded 
regions represent 67% and 90% confidence limits of /. Also plotted are the 790 for which / = 0.5 at each 
data set, the numerical values are given in table 1 (Bromberg et al., 201 Id). 
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We restrict the analysis of Swift and Fermi data to T90 > 0.06 and 0.3 s respectively due to lack of statistics. 



6. Discussion 

The observed plateaus in all three duration distributions, and most notably in the distribution 
of soft BATSE bursts, provide a direct support for the Collapsars model for LGRBs. An inspection 
of different regions of the observed temporal distribution (Fig. 1), under the interpretation of the 
plateau as an imprint of the time it takes the jet to break out of the envelope, provides further 
important information. 

1. The end of the plateau and the decrease in the number of GRBs at long durations, allow us 
to estimate the typical time it takes a jet to breakout of the progenitor's envelope. All three 
distributions are flat below ~ 10 s in the GRB frame, implying that % ~ a few dozen seconds. 



This fits nicely with the canonical GRB parameters used in Eq. |2.1| , and provides another 
indication that the stellar progenitors of Collapsars must be compact (Matzner, 2003). 

2. The distribution at long durations can be used to set a lower limit on the number of chocked jets, 
by extrapolating the slope at Tgo » % to durations shorter than % (see section [?]). At durations 
?90 ^ 100 s Py can be fitted by a power law, p r oc with — 4 < a < —3. Extrapolating 
to Tgo < 4 we find that if p e continues with this power law to t e < tt there are significantly 
more chocked GRBs than long ones. For example, even if we extrapolate this distribution only 
down to t e = 4/2 there are ~ 10 times more chocked GRBs than long GRBs. This prediction 
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is consistent with the suggestion that shock breakout from these choked GRBs produces low 
luminosity smooth and soft GRBs (Nakar & Sari, 2011, and references there in). Indeed the 
rate of such low luminosity GRBs is far larger than that of regular long GRBs (Soderberg et al., 
2006). 

3. At short durations, the duration distribution is dominated by non-Collapsars SGRBs. These 
bursts are hard to classify and are commonly defined by their shorter duration. Our analysis 
shows that putting the dividing line between Collapsars and non-Collapsars at 2 s is statistically 
reasonable for BATSE, and possibly also for Fermi bursts. However, it is clearly wrong to do 
so for Swift. We also calculate the probability as a function of T^q that a burst with a given Tgo 
is not produced by a Collapsar. 

4. While the difference in the lower limit of the plateaus is understood qualitatively as it depends 
on the spectral range of the detector, the variance in the upper limit is less obvious. It may 
reflect various selection effects in triggering algorithms. A more interesting possibility is that 
it reflects a physical origin, e.g. different satellites probing populations with different %. This 
could be explored when the statistical sample of Fermi GBM becomes sufficiently large. 
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